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We present the design of a structured material supporting complete absorption of sound with a broadband 
response and functional for any direction of incident radiation. The structure which is fabricated out of 
porous lamellas is arranged into a low-density crystal and backed by a reflecting support. Experimental 
measurements show that strong all- angle sound absorption with almost zero reflectance takes place for a 
frequency range exceeding two octaves. We demonstrate that lowering the crystal filling fraction increases 
the wave interaction time and is responsible for the enhancement of intrinsic material dissipation, making 
the system more absorptive with less material. 



The damping of sound waves can be understood by conversion of the mechanical energy into heat. The 
acoustic equivalent to an ideal black body would be something similar to a "deaf body which is an object 
absorbing sound coming from all directions at any given frequency. Absorption of sound waves is governed 
by the effects of viscosity and thermal conduction in fluids. In order to give a macroscopic picture of these 
dissipative processes one must interchange the mass density p and the bulk modulus K with complex quantities, 
leading to a finite penetration length into the dissipative medium. In a first approximation, damping can be 
neglected for audible sound in free space and therefore the utilization of lossy materials is essential for screening 
noise and designing acoustic insulators. Thus, it is desirable to create structures with the capability of efficient 
dissipation, preferably in a way that energy conversion causes all waves to be absorbed such that no back- and 
through-radiation takes place. While this remains to be the ideal case, various composites and artificial structures 
were designed in the attempt of pursuing this ultimate goal. We distinguish between resonance-based and 
broadband systems. Locally resonating materials have been fabricated in the form of mass-loaded thin mem- 
branes, gas-bubble arrays and elastic beams, featuring sharp and narrow absorption peaks 1 " 5 . Broadband absorp- 
tion, on the other hand, has been demonstrated for low frequencies by lattices of perforated shells 6 . Sound 
blocking screens made out of, e.g., phononic crystal or metamaterials with one single negative effective parameter 
prevent waves to penetrate through these structures resulting in full reflection of sound 7 " 10 . These effects are 
caused by diffraction properties and evanescent modes, respectively, and will not easily lead to broadband 
absorption of sound. Periodic penetrable structures, generally speaking, have been fabricated with many facets 
for different kinds of waves. Electromagnetic (EM) structured materials, in that regard, such as gratings with finite 
conductance, convex grooves or the moth eyes, comprise anti- reflective systems with a broad spectral res- 
ponse 11 " 15 . However, to sustain a spectrally and angularly rich performance with complete absorption and little 
material use remains a challenge worth pursuing. 

We present a system inspired by recent EM experiments where a forest of vertically aligned single-walled 
carbon nanotubes showed extremely low reflectance and making it the darkest man-made material ever 16 20 . The 
mechanism consists in the attempt of matching the material index to free-space to prevent back- reflected waves 
but at the same time providing sufficient material losses to guarantee intensity attenuation. Here, we show how a 
low- density porous lamella array, in analogy to its electromagnetic counterpart made of nanotube arrays, behaves 
most closely like a true deaf body. Within this framework, we show how these constructed crystals in a most 
counter intuitive way become more absorptive when less material is chosen. 
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Figure 1 | Controlling the absorption by lowering the crystal filling 
fraction f. (a), Schematic of the crystal made of lamellas of width W, lattice 
constant a and length L. The crystal is backed by a rigid support into which 
no sound waves penetrate. (b,c), The enhancement factor 7 and the 
absorption are plotted for various filling fractions with constant lamella 
width W = 4 cm and length L = 0.5 m. In all cases, we implemented 
complex dispersive material dependence for p and K within the lamellas 22 . 

Results 

Concept of enhanced absorption in lossy lamella-crystals. To 

illustrate this surprising behaviour, we begin by analysing a simple 
asymmetric system consisting of an inhomogeneous lossy medium of 
length L supported by a rigid backing as illustrated in Fig. la. The 
general complex scattering matrix can be written as: 

~r(co) 01 
0 1 



S(co) 



(1) 



where r(w) is the complex reflection coefficient. In the absence of 
losses the scattering matrix is unitary and consequently \r(co)\ = 1, 
i.e. perfect reflection due to the rigid backing. This changes 
dramatically in the presence of weak absorption provided that the 
time delay is comparable to 1/21". To see this explicitly, we follow 
considerations developed in the context of absorption in chaotic 
cavities 21 . For a weak absorption rate T the full scattering matrix S 
is related to its loss-less counterpart S 0 



S(w) 



-So(co)(1--Q(cd] 



N 

ffi 4- 

o 

CD 



0.25 0.50 0.75 1.00 
Length L (m) 



% 
o 

in 
43 



0.95 





1 


— 1.3 kHz 

— 2.0 kHz 

— 3.8 kHz 

— 4.0 kHz 


/ 




n 




Figure 2 | Acoustic response of crystals made out of porous lamellas. 

(a), Absorption versus frequency and lamella length L for normal incident 
radiation where f = 30% and W = 4 cm. (b,c) Simulated angular 
dependence of the structure in (a) with L = 0.5 m within the first Brillouin 
zone. We take illustrative examples corresponding to four frequency cuts 
in (c) and show the dependence of absorption with parallel momentum k x . 



where Q(co) = — iSl(co) 



dSpjcp) 
dco 



is the Wigner-Smith delay-time 



matrix. Since S 0 is Hermitian so is Q and its positive eigenvalues 
are interpreted as the group delay time t for a wavepacket centered 
around co. Conservation of energy in the asymmetric system requires 
|r(co)| 2 + A = 1, where A is the absorption, hence, to a first 
approximation we can write 



| r (co)| 2 = l-2rT. 



(3) 



(2) 



Provided a sufficiently long time delay, the reflectance can be 
minimal and the absorption perfect, even for a modest intrinsic 
absorption rate. The perspective of acoustic functional materials is 
to support delay times much exceeding the delay time of their bulk 
counterparts (t^>t/). For a bulk material lamella (/) of length L, 
backed by a perfect reflector, the acoustical path becomes 2L and 
consequently the time delay is X\ = IL/q where q is the sound speed 
in the lamella material. On the other hand, composing a porous 
periodic crystal from the same base material, q is replaced by the 
actual group velocity v g associated with the phononic band 
dispersion relation. In this context, we define the enhancement 
Ft 

factor 7= — — that expresses the acoustic interaction strength in 
the crystal. In other words, enhancement is reached whenever 
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Figure 3 | Measured and simulated absorption and reflectance demonstrating increased dissipation when lowering the filling fraction, (a), Bottom 
view on the suspended crystal made out of porous foam lamellas and supported by a rigid backing. The lamellas have the width W = 4 cm and length L = 
0.5 m. (b,c), Measured (open circles) and simulated (full lines) absorption and reflectance versus frequency for three different filling fractions, f = 40%, 
70% and 100%. All data are obtained for a normal incident acoustic plane wave. 



7 > 1, meaning that dissipation inside the crystal exceeds the 
intrinsic material losses of the lamellas. From Eq. (3) it follows that 
reflectance (absorption) will decrease (increase) whenever wave 
slowing or increased dissipation take place. Since we will not 
operate within the long wavelength regime, we apply this model in 
our interpretation of the enhanced absorption because homogeni- 
zation of the lamella-crystal cannot be applied safely. 

We examine the structure depicted in Fig. la where the goal is to 
explore the tunability for maximal absorption. The basic structure is 
a ID sonic crystal consisting of an array of lossy slices supported by a 
perfectly rigid backing. With translational invariance in the y dir- 
ection we numerically simulate the complex wave interaction by 
coupling free-space sound radiation to Bloch states inside the crystal 
(detailed derivations are found in the Methods section). This leads to 
the following expression for the reflection coefficient 

r G = 2 pfBjeos qj L - S 0G , (4) 
j 

where p ; G , Bj and qj are the/th eigenvector, modal amplitude and out- 
of-plane wave vector of the crystal respectively, from which we derive 
the overall absorption 



where k^ = y kl~ (k x + G x ) 2 . We begin by analysing the spectral 
absorption A (Eq. (5)) for sound incident along the out-of-plane 
crystal orientation at normal incidence. In the experimental set-up 
we utilize a porous material based on homogeneous foam, hence to 
conduct numerical simulations we have implemented dispersive 
material dependence for these lossy materials and calculated the 
acoustical response of the crystal to an incoming sound wave 22 . In 



order to lower the filling fraction and by this controlling the overall 
dissipation, we need to change the lattice constant for a fixed lamella 
width. By inspecting the bulk properties for the case when f = 100% 
or equivalently 7 — 1 we predict overall weak absorption for low 
frequencies as seen in Fig. lc. Obviously this is the regime to be 
challenged since porous materials naturally perform most efficient 
at higher frequencies. Changing the crystal filling fraction f by vary- 
ing the lattice parameter a causes the lattice singularities to shift 
accordingly as seen in Fig. lc. These spectral dips resulting in low 
absorptions arise from diffracted waves becoming grazing and take 
place when \k x + G x \ = co/c 0 23 . If the filling fraction is lowered, a 
growth in the enhancement factor is seen, which is giving rise to 
increased absorption. This increase in dissipation is simply explained 
by the fact that sound is trapped more efficiently inside the crystal as 
compared to the bulk (f = 100%) meaning that 7 must be larger than 
unity as depicted in Fig. lb. We obtain enhanced absorption caused 
by an increased acoustic interaction strength (7 > 1) that for higher 
frequencies slightly grows inversely with the filling fraction, see 
Fig. lb. On the other hand, one could in as much create systems with 
complete absorption 7 > 1 of sound by designing lamellas made out 
of an intrinsic lossier material. It is the choice of materials with high 
absorption rates and the resulting delay time within the dispersive 
crystal which dictates the performance. 

Further to this, we simulate the dependence of the absorption with 
lamella length L. At a given frequency the waves penetrate into the 
crystal with a characteristic length that should not surpass the crystal 
length L > l p to avoid strong back- reflection at the rigid support. For 
this reason, either the wave needs to decay rapidly or travel a suf- 
ficient long distance to undergo enough intensity attenuation as 
plotted in Fig. 2a, showing that absorption increases when the struc- 
ture becomes lengthy. In Fig. 2a we also observe some fine oscilla- 
tions in the mapping of which the spectral locations scale with 1/L as 
known from cavity resonances. More importantly however, although 
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we predict high absorptions in the range 0.95-0.999 in a representa- 
tive band between 2.0 to 2.5 kHz for L = 0.1 m, as seen in Fig. 2a the 
bandwidth of complete sound absorption is easily broadened by 
increasing the lamella length L. Unlike membrane type absorbers 
decorated with rigid platelets where the absorption is resonance- 
based and overall narrow in width 1 , the present scheme relying on 
diffraction constitutes an acoustically thick layer which is essential 
for broadband applications. Next, we investigate the angular sens- 
itivity for the crystal, which up to this point has been simulated for 
normal incident acoustic plane waves only. For this, it is useful to 
calculate the absorption versus frequency and parallel momentum k x 
within the first Brillouin zone. Fig. 2c shows this band diagram con- 
taining Bragg-folded sound lines which are responsible for the lattice 
singularities of low absorption. The entire angular and spectral res- 
ponse resembles the one seen with bandgap materials with no mater- 
ial contrast, but for the present case comprising out-of-plane 
propagation, it is shown that near to 100% absorption sustains for 
the angles shown within the first Brillouin zone (see Fig. 2b for 
specific frequency cuts). Beyond this zone, absorption remains high 
up to grazing incidence and the limits are met only whenever irra- 
diation is in phase with the crystal lattice and full reflection takes 
place. 

Experimental verifications. The crystal is constructed out of thin 
layers of homogeneous foam that is created by compaction and 
compression of a polyurethane, polyester and polyether mixture. 
These lamellas have a width of W = 4 cm and are arranged into a 
ID lattice, which is suspended and backed by a rigid wood support as 
illustrated in Fig. 3a. We have applied the transfer function method 
consisting of a loudspeaker and two microphones for the phase and 
amplitude measurements to detect complex reflections. In the 
experimental set-up the loudspeaker is placed at a sufficient 
distance from the microphones and the sample to ensure plane 
wave generation for all relevant frequencies (see Methods and 
supplementary material). Spectrally, we evaluate the reflectance 
|r(co)| 2 and the resultant absorption A over frequency ranges 
relevant to road and air traffic-noise screening. From the previous 
study we predicted that lowering the crystal filling fraction will 
improve the absorption of sound due to the enhancement of the 
interaction strength. To validate this experimentally we have 
constructed various samples made out of the same lamellas but 
varied the size of the unit cell resulting in three different filling 
fraction as seen in Fig. 3. Bulk material properties are obtained 
measuring the response of the structure with a filling fraction of 
100% from which relatively strong absorption stems from intrinsic 
material losses and the slab length L. When lowering the filling 
fraction down to 70% and 40%, we observe increased perfor- 
mances with mean absorption of 0.97 and 0.99, respectively, over 
an extended spectral range spanning from 0.7 to 3 kHz. Due to a 
finite number of unit cells (see Methods) we detect oscillations in the 
spectrum, overall however, the theory agrees very well with the 
average absorption evaluated from experiments. The crystal 
fabricated with a filling fraction of 40% does not only outperform 
the other samples measured, see insets of Fig. 3, but further to this we 
measure a reflectance on the order of 10" 4 in a representative 
frequency range. 

Due to the interaction strength exceeding unity, complete sound 
absorption extends much further away from normal incidence for 
various directions (k x = k o sin0) as rendered in Fig. 2c. We conducted 
angular- resolved absorption measurements where we vary the angle 
0 from the normal to the crystal surface. Caused by diffraction assoc- 
iated by momentum transfer to the lattice (dash-dotted lines), again 
we predict regions of higher reflections as seen in Fig. 4. This narrow 
region is however being surpassed by a spectrally broad region of 
strong absorption spanning from 6 = 0°-40°. Both, the finite num- 
ber of unit cells and the rigid rotatable frame used for the angular 
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Figure 4 | Measured and simulated angle-resolved absorption spectrum. 

(a), Measured (see Methods for details) and (b), simulated absorption as a 
function of frequency and angle of incidence for the same structure as in 
the previous measurements, now just with f = 36%. The dash-dotted line 
on both contour-maps indicate the condition where I k x + G x \ = w/c 0 . 
(c), Angular dependence of the absorption, representing three measured 
(open circles) and simulated (full lines) frequency cuts from the contours 
in (a) and (b). 



measurements (see Methods) are causing additional unwanted 
reflections in the measurements, Fig. 4a. However, upon inspecting 
both experimentally and numerically the absorption for three differ- 
ent frequencies within the entire resolved spectra (Fig. 4a and 4b, 
respectively), we find overall good agreement validating broadband 
absorption for almost any direction as seen in Fig. 4c. 

Discussion 

We have shown a new concept to engineer highly absorptive acoustic 
materials by means of increasing the sound-material interaction 
strength which is leading to vanishing reflectance, hence, perfect 
absorption of sound. We have achieved this goal by fabricating crys- 
tals made out of intrinsically lossy lamellas and mounted them onto a 
rigid backing. For a representative broad range of frequencies we 
show that absorption can be controlled by the crystal density and 
tuned up to 99%. The counterintuitive ability to increase absorption 
by lowering the effective amount of material is providing many 
interesting opportunities for minimizing noise pollutions by pro- 
ducing efficient acoustic sealing. 
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Methods 

Numerical simulation. We have implemented a plane wave expansion technique 
comprising in-plane Bloch waves and out-of-plane harmonic pressure waves 
p = e lcot e iqz ^2p G e '( kx + G ^ x where q is the out-of-plane wave vector, k x + G with 

G 

2n A 

G = G x = n — x, the in-plane wave vector and p G the eigenvector. After expanding the 
a 

fields and the complex material parameters p and K we end up with the following set 
of equations: 



E 



1 



Pg-c 



(k x + G)-(k x + G') + q 



1 



Pg-o 



Kg-g 



(6) 



Since external radiation impinges onto the crystal that is backed by a rigid support 
(see Fig. la), we need to treat those two regions separately. We write the free-space 
wave radiation in the form of an incident and back- scattered wave 



(k x + G)x-ik G z z 



and the crystal cavity modes as 

p" = 2 E Y.pj B i cos <?;(*-iy (fc+G,x 



(7) 



(8) 



In this cavity, the wave amplitudes are weighted over the j'th eigenvector pj and 
modal amplitude Bj as shown in Eq. (8). When imposing continuity of the fields at the 
interface, we solve the unknows as for example the reflection coefficient r G and the 
overall absorption A (Eq. (4) and Eq. (5) respectively). 

Experimental set-up. The crystal, which is backed by a rigid wood panel, contains 9 
unit cells of lamellas of width W = 4 cm, length L = 50 cm and depth D = 100 cm 
that are mounted upside down without weight suspensions, see Fig. 3a. For the angle- 
resolved measurements, the crystal was mounted on a rotatable frame with a vertical 
axis. In order to fit the 9 unit cells inside this frame, the crystal was constructed with a 
filling fraction off= 36%. 

As rendered in Fig. SI (supplementary material), the experimental set-up con- 
structed to measure the complex reflection coefficient r(co) is based on the transfer 
function method and applied in the anechoic chamber. The first microphone has been 
placed at a sufficiently large distance to the source to ensure propagation of plane 
waves for the entire range of frequencies analysed. Moreover, to avoid finite size 
effects of the truncated crystal, the second microphone is placed in the nearest 
proximity to the crystal. For the sound excitation we use a Genelec 8130A Digital 
Monitoring System placed with a distance of d =4 m to the first microphone. Two 
B&K 4189 microphones of 1/2 inch are placed at the equatorial plane of the crystal 
and aligned with the central unit cell. Data from the two microphones are acquired 
with a CLIO 10 FW (Audiomatica) and then processed in a computer. 
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